Melting of the Electron Glass 
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A model of spinless interacting electrons in presence of disorder is examined using an extended 
dynamical mean-field formulation. When the interaction strength is large as compared to the Fermi 
energy, a low temperature glassy phase is identified, which in our formulation corresponds to a 
replica-symmetry breaking instability. The glassy phase is characterized by a pseudo-gap in the 
single particle density of states, reminiscent of the Coulomb gap of Efros and Shklovskii. Due to 
ergodicity breaking, the "zero-field cooled" compressibility of this electron glass vanishes at T=0, 
consistent with absence of screening. When the Fermi energy exceeds a critical value, the glassy 
phase is suppressed, and normal metallic behavior is recovered. 
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Understanding the basic physical processes determin- 
ing whether a material is a conductor or an insula- 
tor continues to be a central theme of Condensed Mat- 
ter Physics. In recent years, much progress has been 
achieved in understanding weakly disordered metals [Q, 
but the physics of the metal-insulator transition (MIT) 
in presence of disorder remains largely an open problem. 

The interplay of the electron-electron interactions and 
disorder is particularly evident deep on the insulating 
side of the MIT. Here, both experimental and theo- 
retical studies Q have demonstrated that they can lead 
to the formation of a soft "Coulomb gap" , a phenomenon 
that is believed to be related to the glassy behavior of the 
electrons. Such glassy freezing has long been suspected 
to be of importance, but very recent work |^ has sug- 
gested that it may even dominate the MIT behavior in 
certain low carrier density systems. The classic work of 
Efros and Shklovskii has clarified some basic aspects of 
this behavior, but a number of key questions remain. In 
particular: (1) What is the precise relation of this glassy 
behavior and the emergence of the Coulomb gap? (2) 
What should be the order parameter for the glass phase? 
(3) How should the glassy freezing affect the compress- 
ibility and the screening of the electron gas? (4) How do 
the quantum fluctuations (electron tunneling) melt this 
glass and influence the approach to the MIT? 

In this Letter, we examine a simple model for which all 
these questions can be rigorously answered in the limit of 
large coordination. Within our model, the universal form 
of the Coulomb gap proves to be a direct consequence of 
glassy freezing. The glass phase is identified through the 
emergence of an extensive number of metastable states, 
which in our formulation is manifested as a replica sym- 
metry breaking instability. As a consequence of this er- 
godicity breaking, the zero-field cooled compressibility is 
found to vanish at T=0, suggesting the absence of screen- 
ing in disordered insulators. Finally, we show that quan- 
tum fluctuations can melt this glass even at T = 0, but 



that the relevant energy scale is set by the electronic mo- 
bility, and is therefore a nontrivial function of disorder. 

In order to concentrate on the physics of the electron 
glass, we focus on a simple lattice model of spinless in- 
teracting electrons in presence of randomness, as given 
by the Hamiltonian 



(1) 



Here, Ui = clci are the occupation number operators, 
tij are the corresponding hopping elements, Vij describes 
the inter-site electronic repulsion, and Si are the random 
site energies. Motivated by the long range nature of the 
Coulomb interaction, we wish to examine the situation 
where a given electron interacts with a large number Zy 
of other electrons. The analysis is simplest if the inter- 
action Vij — V is taken to be uniform within a volume 
containing Zy neighbors, as opposed to the more realistic 
Coulomb interactions. Nevertheless, most of the qualita- 
tive features of the Coulomb glass are still captured if we 
consider Zy >> 1, as follows. 

To obtain the relevant dynamical mean-field (DMF) 
equations ||^, we focus on a given lattice site, and for- 
mally integrate out |Q all the other degrees of freedom. 
For large coordination, the corresponding local effective 
action simplifies, since it suffices to consider contribu- 
tions up to the second order in the hopping elements t 
and the interaction amplitudes V We concentrate 
on a Bethe lattice at half filling, where the equations are 
simplest, and we find 
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We have used standard functional integration over 
replicated Grassmann fields where a — l,...n 

(n 0) are the replica indices Here, the oper- 

ators Snf{T) = (cJ"(t)c°(t) — 1/2) represent the den- 
sity fluctuations from half filling. The order parameters 
G{t — r'), x{t — t') and qab satisfy the following set of 
self-consistency conditions 

G(r - t') = / de,P{e,) < cl''(T)cnr') >eff, (3) 

X(r -r')=J de,P{e,) < <5nl"(r)K(T') >eff, (4) 

qab = / d£^P(£^) < fcj'^ (r)<5n,^ (t') >e// . (5) 



In these equations, the averages are taken with respect 
to the effective action of Eq. (2), and P{ei) is the prob- 
ability distribution of random site energies. Physically, 
the quantity G(r — r') represents the single-particle elec- 
tronic spectrum of the environment, as seen by an elec- 
tron on site i. The second quantity x(t — t') is the aver- 
aged dynamic compressibility. This mode- coupling term 
reflects the retarded response of the density fluctuations 
of the environment. Finally, the quantity qab {a 7^ b) is 
the familiar Edwards- Anderson order parameter pl . Its 
nonzero value indicates that the time averaged electronic 
density is spatially non-uniform. 

We should emphasize that in presence of randomness, 
the electronic density will remain nonuniform at any tem- 
perature. As a result, the Edwards- Anderson order pa- 
rameter q = < Sn ^ everywhere, and thus cannot 
be used to identify the glass transition. However, in the 
classical (t — 0) limit, a simple transformation Si — 2dni, 
J = V/A, and hi — ei/2 reduces our problem to the fa- 
miliar random-fleld Ising model. Here, recent work jT^ ] 
has suggested that a glassy phase can be identified at 
low temperatures by carrying out an appropriate replica 
symmetry breaking (RSB) stability analysis. 

To carry such a RSB analysis in our approach, note 
that the our self-consistency conditions Eqs. (2)-(5) can 
be derived from an equivalent infinite range model given 
by Eq. (1), but having random hopping elements tij and 
random inter-site interactions Vij [ p^ , with zero mean 
and variance t'^/N, and V'^/N, respectively (iV is the 
number of lattice sites). This model is similar to familiar 
spin-glass models, and the RSB analysis can be carried 
out using standard methods The resulting stability 
criterion takes the form 



d£iP{£i) [xioc{£i)Y 



(6) 



Here, xioc{£i) is the unaveragedlocaX compressibility, that 
can be expressed as 



d 1 



/3 



X^oc(£^) ^ g^^ I d-T < Sn{T) >Se//(») ■ (7) 



We continue by concentrating on the classical [t = 0) 
limit, where x(t, r') = 1/4, and the problem can easily 
be solved in closed form. We first focus on the replica 
symmetric (RS) solution, which is valid at high temper- 
atures, and set qab = q for all replica pairs, giving 



1 



q = 



+ 00 



dx 



e--'/2tanh2 



xPWeffiq) 
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where we have introduced an effective disorder strength 
Weff — \JW^ -I- V^q , and have considered a Gaussian 
distribution of random site energies of variance . Even 
for weak disorder, g 7^ at any temperature, but it 
becomes appreciable only for T < V/4. These frozen- 
in density fluctuations introduce an added component 
to the random potential seen by the electron. In this 
way Weff can become appreciable at low temperatures 
even if the "bare" disorder W is weak; this represents 
an interaction-driven mechanism for resistivity enhance- 
ment. 

Next, we examine the stability of this RS solution. 
From Eq. (6) we flnd an expression for the instability 
line, which takes the form 
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, (9) 



where q{T,W) is given by Eq. (8). The resulting ex- 
pression for the glass transition temperature simplifles 
for W >> V, where using the fact that q < 1/4 we 
flnd Tg « ^=w- This reduction of To at large W 

Dv27r 

reflects the reduction of frustration in the strong disor- 
der limit, where the electrons simply populate the lowest 
potential wells. Still, in contrast to the well known "AT- 
line" ||ll|, Tg decreases slowly with disorder, suggesting 
that the glassy behavior of electrons may be observable 
at finite temperatures, in agreement with some recent 
experiments |l4[ . 

To understand this behavior, we investigate the struc- 
ture of the low-temperature glass phase. Consider the 
single-particle density of states at T=0, which in the clas- 
sical limit can be expressed as 



(10) 



where ef' = Si -\- ^ij^j are the renormalized site ener- 
gies. In the thermodynamic limit, this quantity is noth- 
ing but the probability distribution PR(ef ). It is anal- 
ogous to the "local field distribution" in the spin-glass 
models, and can be easily shown to reduce to a sim- 
ple Gaussian distribution in the RS theory, establishing 
the absence of any gap for T > Tg- Obtaining explicit 
results from a replica calculation in the glass phase is 
more difficult, but useful insight can be achieved by us- 
ing standard simulation methods 1^,^ on our equiva- 
lent infinite-range model; some typical results are shown 
in Fig. 1. 
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FIG. 1. Single particle density of states in the classical 
{t — 0) limit at T = 0, as a function of disorder strength. 
Results are shown from a simulation on A*' = 200 site sys- 
tem, for W/V = 0.5 (thin line) and W/V = 1.0 (full line). 
Note that the low energy form of the gap takes a universal 
form, independent of the disorder strength W. The dashed 
line follows Eq. (11). 

We find that as a result of glassy freezing, a pseudo-gap 
emerges in the single-particle density of states, reminis- 
cent of the Coulomb gap of Efros and Shklovskii (ES) 
The low energy form of this gap appears universal, 



p{e)^C\er/V'; C 



1, 



(11) 



independent of the disorder strength W, again in strik- 
ing analogy with the predictions of ES ||l^ . To establish 
this result, we have used stability arguments very sim- 
ilar to those developed for spin-glass (SG) models [0, 
demonstrating that the form of Eq. (11) represents an 
exact upper bound for p{e). For infinite-ranged SG mod- 
els, as in our case, this bound appears to be saturated, 
leading to universal behavior. Such universality is often 
associated with a critical, self-organized state of the sys- 
tem. Recent work [ pT| finds strong numerical evidence 
of such criticality for SG models; we believe that the 
universal gap form in our case has the same origin. Fur- 
thermore, assuming that the universal form of Eq. (11) 
is obeyed immediately allows for an estimate of Tg{W). 
Using Eq. (11) to estimate the gap size for large disorder 
gives Tg ^ Eg ^ V'^ /W, in agreement with Eq. (9). 

The ergodicity breaking associated with the glassy 
freezing has important consequences for our model. 
Again, using the close similarity of our classical infinite 
range model to standard SG models [0, it is not diffi- 
cult to see that the zero- field cooled (ZFC) compressibil- 
ity vanishes at T = 0, in contrast to the field-cooled one, 
which remains finite. Essentially, if the chemical poten- 
tial is modified after the system is cooled to T = 0, the 
system immediately falls out of equilibrium and displays 
hysteretic behavior [p^ with vanishing typical compress- 
ibility. If this behavior persists in finite dimensions and 



for realistic Coulomb interactions, it could explain the 
absence of screening in disordered insulators. We em- 
phasize that our picture does not agree with the recently 
proposed scenario of Si and Varma [T^ , where the two 
compressibilities are not distinguished, and are both pro- 
posed to vanish in the T — Q insulating state. 

Finally, we examine the effect of quantum fluctuations 
on the stability of the glass phase. We generally expect 
that for a sufficiently large Fermi energy Ep ^ to, Tq 
0, but of particular interest is the disorder dependence of 
tciW). In the quantum case, solving our DMF equations 
becomes a non-trivial many-body problem, which cannot 
be solved in closed form. The problem simplifies in the 
strong disorder limit, where the exact expression for the 
RSB boundary can be obtained. As in the classical case, 
to examine the W >> V limit, we use the fact that 
<Z < 1/4, x{t — '''') ^ 1/4, so that in solving Eq. (6) we 
can eliminate the corresponding terms in the expression 
for the local effective action Eq. (2). The calculation 
of the local compressibilities of Eq. (3) then reduces to 
computing the same quantity for noninteracting electrons 
in presence of disorder, and we find 

Xioc{e) = ^J dLup^{Lu)cosh-'{-PLu). (12) 

Here, the local density of states (LDOS) reduces to a 
narrow resonance of width A = TTt'^P{ei = 0) '--^ t^ jW 
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TT (w - e)2 -f A2 ' 



and we get 



taiT = 0, c5o) = V/TrV2. 



(13) 



(14) 



In contrast to the effect of finite temperatures, wc find 
that a finite value of the Fermi energy Ep ^ t is needed 
to suppress the glass phase, even ior W >> V I This 
surprising result has in fact a simple physical interpre- 
tation. At i = we have argued that the glass phase 
for T4^ >> V is characterized by a small energy scale 
Eg ~ V'^ /W, hence the reduction of Tq at strong disor- 
der. In the T — quantum case, the role of the thermal 
energy is replaced by an energy scale measuring the size 
of quantum fluctuations. However, the relevant energy 
scale is not the Fermi energy Ep ~ t, but instead is given 
by the hybridization energy A ~ t'^/W which, according 
to Fermi's golden rule, determines the inverse lifetime of 
an electron on a given site. Due to the decrease of the 
electronic mobility at large disorder, A (14/^) is reduced at 
W >> t; comparing it to the T = gap size Eg{W) 
immediately explains Eq. (14). 

Our solution, which evaluates the local compressibili- 
ties in Eq. (6) by setting V — gives the exact behavior 
for W >> V. This approximation can be considered 
as a leading term in a "weak-coupling" approach, where 
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the compressibilities are expanded in powers of the inter- 
action, similar in spirit to the RPA of the electron gas. 
For our purposes of calculating the phase diagram, this 
strategy will provide qualitatively correct results even 
for weaker disorder, and the resulting phase diagram is 
shown in Fig. 2. Note that our glass phase persists even 
as W — *■ 0. However, in this regime we expect uniform 
charge ordering to preempt any glass formation, an ef- 
fect that we have not examined, but that can easily be 
incorporated in our framework. 



freezing strongly enhances the effective disorder seen by 
the electrons, providing a new driving force for electron 
localization. In addition, the ergodicity breaking in the 
glassy state leads to the universality of the Coulomb gap, 
as well as the incompressible behavior in disordered in- 
sulators. 
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FIG. 2. Phase diagram as a function of temperature T, 
disorder strength W, and hopping element t. Note how the 
glass transition temperature To decreases as 1/W in the 
strong disorder limit. In contrast, the critical value of the 
hopping element to remains finite as W ^ cxd. 

Our results clearly demonstrate that the quantum fluc- 
tuations become a nontrivial function of disorder in this 
electronic context, which plays a key role in controlling 
the glassy phase. We should emphasize that, within our 
theory, the relevant hybridization energy A remains finite 
at any t ^ 0, since Anderson localization is absent 
at Zt = oo. As a result, while the glass phase will per- 
sist for all t < t(3, the systems will become metallic for 
any t ^ 0. Similarly as finite temperatures, A ^ will 
introduce a cutoff for both our Coulomb gap, and the 
corresponding ZFC compressibility. On the other hand, 
Anderson has demonstrated |l^] that A in an An- 
derson insulator, even in absence of interactions. We thus 
expect that Anderson localization will strongly stabilize 
the glass phase, as well as restore the ZFC incompress- 
ibility and the existence of the Coulomb gap in the entire 
insulating phase. Our theory can be readily extended to 
incorporate the Anderson localization effects within the 
DMF framework, following the approaches of Ref. [ pO| , 
but these interesting issues will be addressed elsewhere. 

To conclude, we have presented a simple theoretical 
framework where several aspects of glassy freezing of elec- 
trons can be investigated in detail. We find that glassy 
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